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Combining a recent derivation of the CMB evolution equations for homogeneous but anisotropic 
(Bianchi) cosmologies with an account of the full linearized dynamical freedoms available in such 
models, I calculate and discuss the various temperature and polarisation anisotropy patterns which 
may be formed. Certain anisotropics can be hidden in superhorizon modes at early times, thus 
avoiding any constraints from nucleosynthesis while nevertheless producing non-trivial redshift-zero 
temperature patterns in flat and open universes. The results are likely to be more of pedagogical 
than observational interest, but future work will assess whether such patterns can be matched to 
anomalies in WMAP results. 
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I. INTRODUCTION 

Bianchi universes - a large and almost complete class of 
relativistic models which are homogeneous but not neces- 
sarily isotropic - have remained of interest to mathemat- 
ical cosmologists since seminal papers uncovering their 
major properties in the 1950s and 60s [e.g. [H El E] |4]; 
for a concise and accessible review see Ref. [5^. How- 
ever for observers they are often considered obsolete, the 
chaotic cosmology programme having been superseded by 
the inflationary paradigm, especially in the light of CMB 
temperature maps from COBE [6j and WMAP [7 which 
are statistically very close to isotropy. On the other 
hand, there are deviations from the standard model pre- 
dictions, consistent between COBE and WMAP, which 
are not easily linked to known systematics (e.g. [3 |9l [10] 
and references therein). Jaffe et al have shown that 
anomalies such as the low quadrupole amplitude, align- 
ment of low-/ modes and large-scale power asymmetry 
can be mimicked by adding (to the standard, scale-free 
perturbations which dominate the CMB signal) Bianchi 
temperature maps derived for a specific subcase of Yllh 
by Hawking and later Barrow [121 [E] (see also Ref. [14] 
in which the temperature patterns are discussed quali- 
tatively by considering the Killing motions in a local or- 
thonormal frame) . Such models are more predictive than 
simple Bianchi I models which can trivially be fine-tuned 
to solve the quadrupole problem, e.g. \TE\ [TH [17]. Al- 
though the cosmological parameters implied by Jaffe et 
al's Vll/i fit are inconsistent with concordance parameters 
[TSl [T9] and constraints from primordial nucleosynthesis 
[20] , there is a possibility that these tensions may be re- 
solved by a more complete dynamical model [181 Ull [22] • 

Recently, in Ref. ^21j (henceforth PC07), a com- 
plete and computationally convenient Boltzmann hier- 
archy framework for calculating temperature and polar- 
isation anisotropics in any nearly Friedmann-Robertson- 
Walker (FRW) Bianchi universe was described. For the 



favoured parameters based on Jaffe et al's analysis of ob- 
served temperature maps, the predicted B-mode polar- 
isation is too strong to be consistent with observational 
upper limits from WMAP three year data [23]. However, 
although the theoretical analysis was general, in numeri- 
cal results PC07 followed previous Bianchi CMB studies 
in employing dynamical solutions for which the shear de- 
cays as cr cx a~^, where a is the FRW scale factor. It is 
as a direct result of this rapid decay that the polarisation 
amplitude is generally very strong relative to the anoma- 
lous temperature anisotropics; thus, like nucleosynthesis 
and density parameter inconsistencies, the polarisation 
conflict may be resolved by considering more general dy- 
namics. 

In this note, I discuss the appearance of the Bianchi 
CMB in the generalised nearly-isotropic regime. (Models 
which are highly anisotropic at any point after recombi- 
nation are ruled out by observations if one assumes the 
Copernican principle [2l]-) I will show that the afore- 
mentioned polarisation and nucleosynthesis conflicts are 
simultaneously, yet naturally, avoided in a certain class 
of models. By fitting these models to CMB anomalies, 
constraints on cosmological parameters can be produced 
and tested against trusted concordance values; however, 
this requires a detailed statistical framework and is post- 
poned to later work. 

The setup and dynamics of the general case are de- 
scribed in ^ and applied to flat and open models in 
^III For pedagogical interest, I give results from mod- 
els with a variety of cosmological parameters which are 
not necessarily consistent with other constraints. I de- 



scribe the results for a closed universe in pVj although 
the anomalous CMB contribution in such cases is always 
a quadrupole and thus less interesting from a phenomeno- 
logical standpoint. Finally, I summarise the work in 3Vj 



II. SET-UP AND DYNAMICS 
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During an epoch of near-isotropy, the dynamics of 
any Bianchi model can be analyzed by a decomposition 
into an isotropic FRW background and particular linear 
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perturbations which break the isotropy while maintain- 
ing homogeneity with respect to three simply-transitive 
Killing vector fields (KVFs). For example, one might 
take an expanding fiat model and perturb the expan- 
sion rates along one or two of the perpendicular coordi- 
nate axes; this canonical example arises as a special case 
(Type I) in the formalism described below. 

The geometrical setup is described by the structure 

constants Cf, which fix the commutators of the three 

'■J 

preferred Killing fields 



(1) 



C% and the Jacobi 



The explicit antisymmetry C^^ - 
identities C^i^JJf-^^ = together substantially restrict the 

possible values for the C^^-; in fact, up to arbitrary scal- 
ings and rotations only four of the initial 27 values are re- 
quired to describe a given space (canonically [1 these are 
known as ni,n2,n3 and ai, where ai should be sharply 
distinguished from the FRW scale factor a). The cate- 
gorization of Bianchi models into distinct types rests on 
classifying the inequivalent sets of these four constants; 
here I will consider only the small subset which admit 
an FRW limit. Succinct but explicit derivations of the 
relations between structure constants and the resulting 
isotropic background cosmologies are available in PC07; 
or for more detail see e.g. Refs. [27| l29l 130]. 

I will first consider the open and fiat models (Bianchi 
type Vllh and its limiting types I, V and VIIq); these 
possess two distinct geometric freedoms. The first corre- 
sponds to a characteristic length over which the Killing 
fields spiral relative to a parallel-propagated frame. In 
the canonical decomposition the physical scale of this 
spiral is fixed by setting ni = 0, n2 = ns = 1 and scal- 
ing the present horizon size via the dimensionless Hubble 
parameter, usually denoted x [121 |3T]. The second free- 
dom corresponds to the isotropic curvature scale and can 
be set by specifying the fiducial FRW parameter in 
terms of which ai = ztx^/^K- This sign ambiguity re- 
fiects a parity freedom; we will adopt the conventions of 
PC07 noting the standard result that, in terms of the 
CMB, parity inversion can be achieved by mirroring all 
maps and inverting the sign of B polarisation. 

As X ^ oo for fixed fix one pushes spiral structure out 
of the horizon and obtains Type V models (this situation 
can be renormalized so that x = l/v^H^, n2 = na = 
and ai = ±1); for ft^ = we have the flat Type VIIq 
specialisation (n2 = na = 1, ai = 0); with x ^ oo 
and CIk = one obtains Type I (renormalisable such 
that X takes an arbitrary finite value with all structure 
constants vanishing), which was described at the start of 
this section. 

Except in type I, the choice of canonical values pres- 
elects a preferred axis, rotations about which leave the 
structure constants invariant; the alignment of the per- 
turbed shear eigenvectors relative to this axis determine 
both the patterns obtained in the CMB and the na- 
ture of the dynamics, but rotations around the axis are 



pure gauge. Thus anisotropic perturbations are divided 
into three classes which transform respectively as scalars 
(henceforth 5-modes), vectors (I'-modes) and rank-2 ten- 
sors (t-modes) about the locally preferred axis; modes in 
these classes must evolve independently in the linear ap- 
proximation [44]. The modes generate in turn m = 0, 1 
and 2 spherical harmonics in the CMB if the fiducial 
^ = direction is taken along the residual symmetry axis 
(as in PC07). The general nearly-FRW Bianchi CMB ap- 
pears as a linear combination of the three resulting pat- 
terns, although here we will study each mode separately. 

A careful analysis of previous works [25j [26 shows that 
the linear amplitude A of any Bianchi mode evolves ac- 
cording to 



2nA' + = 



(2) 



where primes denote the derivative with respect to con- 
formal time 77, H = a' /a is the conformal Hubble param- 
eter, A is a time-independent complex eigenvalue associ- 
ated with each decoupled mode and the shear resulting 
from the mode has magnitude a oc /a. (When A ^ M 
the argument of the complex-valued A defines the ori- 
entation in the preferred plane; as previously stated the 
evolution equations are invariant under rotations in this 
plane, A Ae^^^^.) Eq. ^ assumes there to be no 
anisotropic matter source terms, which in the absence of 
exotic effects should be a good approximation for the re- 
combination era onwards (even if neutrino viscosity has 
an important role to play during radiation domination 
[e.g. 22 , 28 ). 

Further technical discussion of the Bianchi linearisa- 
tion and its relation to the existing literature will be given 
in a future work [27]; however, its essential features can 
be understood without the extensive framework required 
for a full derivation. For instance one may show that 
each mode maps onto a specific first order FRW metric 
perturbation (e.g. Bardeen [40 ) with A proportional to 
the mode amplitude in the natural synchronous gauge 
defined by the Bianchi-homogeneous space slices. Seen 
from this perspective, the decaying modes correspond to 
vector perturbations (eq. 4.12 - 4.13 of Ref. [40 , with 
Bardeen's ^ = A') and infinite wavelength gravitational 
waves (eq. 4.14 of Ref. [40 , with k ^ and Bardeen's 
H = A) whereas the oscillating modes correspond to 
gravitational waves of specific finite wavelengths (with 
our A = /c^ + 2K^ where K is the FRW curvature param- 
eter and k is the wavenumber in Ref. [40]). The particu- 
lar wavelengths are associated, in the open and fiat cases, 
with the aforementioned spiral lengthscale. Alternatively 
one may discern the physical status of the modes by con- 
sidering gauge invariant covariant variables [30 , in terms 
of which the oscilliatory modes arise from the shear prop- 
agation equation ([30], eq. 31) as an interaction between 
shear and anisotropic 3-curvature {S ex \XA\/a'^). 

For closed (Type IX) models the canonical structure 
constants read ni = 77,2 = na = 1, ai = 0; in these units, 
the Hubble parameter is fixed by the FRW curvature 
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FIG. 1: The v (vector) modes of type VII/i and its specialisations. For a grid of parameters the temperature pattern is shown 
(upper panel) and E and B mode polarisation exaggerated in scale by a factor 20 (small lower left and right panels respectively). 
X varies from top to bottom while Qm (and hence Qk) varies from left to right. I assume Qa = except in the rightmost 
panels for which the concordance values (Qm,^^a = 0.24,0.76) are used. The primary effect of non-zero A is to increase the 
conformal time to last scattering, causing the final z — ^ pattern to be more tightly wound. When Qm + = 1 one obtains 
type VIIo; as x ^ oo one obtains type V (from VII^) or I (from VIIo). The model with (Viyi^x) — (0.5,0.6) (shaded panel) is 
able to mimic known CMB temperature anomalies [111 119] . 



such that X = l/v"— 411k, i-e. = 1/4. If the physical 
curvature radius is made progressively larger, we must 
take VLk and x ^ oo. As before this limit yields, 
as expected. Type I models (explicitly by renormalizing 
such that ni = n2 = na = and x takes an arbitrary 
finite value). Further discussion of the closed models is 
deferred to ffVl 



In all cases the physical value of i^o enters only through 
the recombination history which is obtained using Rec- 
FAST [32 assuming i^o 72kms"^Mpc"^ and 
().() 1 1. but allowing (^m and (^a to vary consistently with 
the Bianchi model. For quantitative results in this peda- 
gogic exploration I have ignored reionisation, discussing 
the effects of such an assumption individually for each 
mode. 



III. RESULTS: OPEN AND FLAT MODELS 



In this section, I will discuss the possible patterns aris- 
ing in the CMB when breaking the isotropy of open and 
flat FRW models; the relevant Bianchi type is VII^ (to- 
gether with its previously described limiting cases). Of 
the three classes of linearized anisotropic freedoms dis- 
cussed above it is the v modes [with A = in Eq. ([2|] 
which have been investigated in previous works on the 
Bianchi CMB [IH [H [H [21]; a grid of such models 
for varying x (vertical) and (horizontal) is shown in 
Fig. [l] A similar plot can be found in Ref. [33 but here 
I include the limiting types (V, VIIq and I) and addition- 
ally display the E and 5-mode polarisation in two small 
panels beneath each temperature map. The model with 
(Vl^.x) = (0.5,0.6) (shaded in Fig. [l]) generates maps 
which mimic known CMB anomalies when suitably ori- 
ented on the sky [TTl [191 EU- The colour scale of the 
temperature plot is arbitrary (since one can rescale the 
linear perturbation to obtain any specified small mag- 
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FIG. 2: The s (scalar) modes of all open and flat types are insensitive to the spiralling of the basis vectors, controlled by x, 
because the quadrupole anisotropy is aligned along the symmetry axis (m = in the decomposition of PC07). Similarly, no 
B-modes can be generated. In the VIIo/I limit, there is very little sensitivity to the value of since the anisotropy is always 
a pure quadrupole; only the relative polarisation strength is marginally affected (because of the increased conformal time to 
last scattering for larger Qa). 



nitude) but the polarisation panels are shown exagger- 
ated relative to the temperature panels by a factor of 20. 
Because the polarisation is already strong, reionisation 
would increase its strength only marginally (see PC07 
for details). 

While most models are calculated assuming A = 0, the 
final column displays concordance (17m, = 0.24,0.76) 
maps. The effect of non-zero A can be understood by 
examining equation (33) of PC07: the final patterns ob- 
tained are determined by ai and the conformal time r] 
to last scattering. Keeping these two quantities constant 
recovers the ((7m,^a,^) degeneracy noted in previous 
works [iHl [T^ (although in the general case this picture is 
complicated somewhat by differing <j(a) behaviour when 
A 7^0). 

The VII/iS mode (Fig. [2]) is similar to the v modes 
in its dynamical behaviour (also having curvature eigen- 
value A = 0) and hence produces similar polarisation am- 
plitudes. Because the mode is invariant under rotations 
about the preferred axis, the familiar spiral structure is 
missing (and hence the value of x has no effect). Fur- 
ther, 5-mode polarisation is not generated (this is nec- 
essarily so when the transformation of the mode under 
reflections leaving the symmetry axis invariant is consid- 
ered). The only remaining effect is the deformation of the 
quadrupole into a focussed spot, which cannot be evaded 
in anisotropic, homogeneous open universes [14[ [27 ] . 

Of more interest are the Wlh^ modes, which are dy- 
namically non- trivial with curvature eigenvalue A = 4(1 — 
iai) in Eq. ([2|. The resulting evolution is obtained nu- 
merically; however analytic approximations are helpful 
in interpreting the results. For instance during matter 
domination, while the isotropic curvature remains on su- 
perhorizon scales, we may assume VLk = and Ti = 2/77, 
yielding the explicit solution 



where the constant k, = a/A/2 = ^/T 
results are unaffected by assuming hz 



2{Kr]f {nr]f 



sin 2tiiri cos 2^7] 
2 {nr]f {nr]f 



- iai but physical 
1 [45;. ^0 and Ai 
are integration constants determined by boundary condi- 
tions, e.g. the shear and anisotropic curvature at a fixed 
time. 



(3) 



Modes with = 0, Ai 7^ {Hi modes') are regular 
as 7^ ^ (i.e. their dimensionless shear a/H = \A'\/H 
and anisotropic curvature S/H'^ = \XA\/H'^ both tend 
to zero); conversely modes with 7^ ('to modes') 
have divergent anisotropy for 77 ^ 0, thus predicting an 
early non-linear phase (which is not itself captured by 
our linear description). Although the picture is some- 
what complicated by neutrino free-streaming, matching 
across matter-radiation equality shows that one may con- 
struct a ti mode with small shear all the way to the initial 
singularity [46 . Thus by choosing a late-time solution 
with Aq = one will automatically satisfy nucleosynthe- 
sis constraints (see also [35 ). Furthermore, if the am- 
plitude of all modes is set by some form of equipartition 
theorem at high redshifts (e.g. [28 ), one may expect on 
entering the linear regime that the divergent to, ^ and 
s modes have similar amplitudes and decay rapidly to 
leave only these regular ti modes. (However numerical 
models show that, with arbitrary initial conditions, the 
amplitude of the ti mode is often sufficient to trigger a 
late-time non-linear anisotropic phase when it enters the 
horizon, e.g. [36 , giving rise to the intermediate isotropi- 
sation picture in which homogeneous shear is only small 
for a finite period Of course if A 7^ 0, this final 

non-linear phase can be avoided by shrinking the hori- 
zon again at late times.) One should bear in mind that 
to modes will generate similar CMB patterns but have 
radically different dynamical behaviour; however, I will 
not consider such a case quantitatively. 

The patterns formed by the ti mode are shown for a 
grid of parameters in Fig. [3j As x ^ 00 the dimension- 
less curvature eigenvalue A/x^ shrinks to zero (and thus 
the anisotropic curvature becomes pure gauge) - so no 
strict Type V or Type I limits exist for ^o = and only 
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FIG. 3: The nucleosynthesis-compatible ti modes of types VII^ and VIIq. The polarisation magnitude is exaggerated in scale 
by the factors specified in each subplot. The primary polarisation varies according to Eq. ^ and, because the shear can grow 
after recombination, is generally much weaker relative to the temperature anisotropy than with s and v modes for which the 
shear always decays. 



finite values of x are considered. As before, the polarisa- 
tion shown is purely primordial and ignores reionisation. 
The large variations in its strength between models can 
be understood by considering the shear near the decou- 
pling surface. For models where ftu is not tiny, effects are 
controlled by the long matter-dominated phase during 
which, unaware of small A or superhorizon isotropic cur- 
vature, observers see a flat Universe with effective Hubble 
parameter ^/Ttux. For the polarisation strength one ap- 
proximates 

Q, (X (TLSS At (X „ , /^rrn 

where At is the expected time between penultimate and 
final scattering, Hq is the physical Hubble parameter and 
/e is the electron fraction at last scattering which scales 
approximately as qI^'^Q^^Hq^ [e.g. |37l; since Hq and 
Q]j have been kept constant for all models, the overall 
polarisation anisotropy is taken to be proportional to 
Aix~'^Q^. On the other hand, the total temperature 
anisotropy is built up between zlss < ^ < 0. Ignor- 
ing the advection of the observed pattern to estimate the 



magnitude of the temperature anisotropy we have 

ATcx / adt= A'dr] = A{r]o)-A{r]Lss) (5) 

where rj^ = 2/ {x^/VLm) is the present day equivalent con- 
formal time. Inserting the solution ([3| one obtains the 
limits 

- ^1 \ 4/3 % » 1. 

The final polarisation-to-temperature ratio, combining 
Eqs Q and (|6|, should scale roughly as 

which agrees with the trends shown in Fig. [3j Unlike 
the s and to modes, here inclusion of reionisation will 
significantly boost the polarisation strength over the pri- 
mordial expectation since the majority of models pos- 
sess a considerably stronger temperature than polarisa- 
tion quadrupole at z < 10. 

An obvious question concerns the existence or other- 
wise of nucleosynthesis-compatible models which, in a 
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similar vein to the original YllhV fits (Fig. [T] shaded 
panel), can mimic known residual CMB anomalies. If 
this were possible using parameters consistent with con- 
cordance values, i.e. with maps from the right-hand col- 
umn of Fig. [3j the type VII Bianchi models would com- 
mand renewed observational interest. This will be inves- 
tigated, comparing the maps developed here with data 
from WMAP in future work. 



IV. COMMENT: CLOSED MODELS 

Type IX models have been shown by Grishchuk, 
Doroshkevich and Yudin [26 to decompose into a closed 
FRW background with superimposed maximal wave- 
length (/c^ = 6K) gravitational wave perturbations (see 
also Ref. [38 for helpful comments) obeying the standard 
propagation equations ([40^, eq. 4.14; see end of pl|). 
Comparison shows that, within our terminology, this is 
modelled by the standard mode evolution Q with A = 2 
(this result is obtained directly from the anisotropic field 
equations in [27]). 

By an extension of the arguments already presented, 
setting K, = V2 in Eq. (|3), one can see that the primary 
polarisation is again small for modes which evade nucle- 
osynthesis constraints; however the temperature patterns 
are less phenomenologically interesting as they remain 
quadrupolar (at first order). 

To see why this is, consider the origin of non-trivial 
patterns in the previously considered open and fiat case. 
The temperature anisotropy AT/T arises from integrat- 
ing the shear tensor cr^j projected, on each timeslice, 
along the tangent vector to the photon path (PC07, 
eq. 22). Taking the normalized spatial components of 
the photon momentum (such that p'^Pi = 1, i = 1 • • • 3), 
one has 

(8) 



oscillates, the microwave background appearing exactly 
isotropic when A{to) = ^(^lss) which may occur an in- 
definite number of times depending on the initial condi- 



where the integral is taken between last scattering (^lss) 
and the time of observation {to) and in general is a 
function both of observation angles Oq^ <po and time t. 

I have implicitly adopted a spatial chart such that all 
fiuid quantities have constant components in a given spa- 
tial slice (this is always possible for Bianchi spaces; see 
PC07). In particular this means that the components of 
the shear tensor a^j depend only on time not on the ob- 
servation angles. Written in this basis, non-quadrupolar 
temperature patterns can only arise through changes in 
the components of p [at zeroth order for a first order so- 
lution to Eq. ([8|]. In the Type IX case the components 
of remain constant to this accuracy [47 ; thus the tem- 
perature anisotropy remains purely quadrupolar. 

According to the definitions in |llj the shear oscil- 
lates in time such that (Jij{t) = a~^' {t)dij for a fixed 
dij. The temperature anisotropy (l8| may be written 
AT/T = (A(to) - A(tLss))^ijPV the fixed obser- 
vation vector p{Oo^(j)o). Thus the quadrupole intensity 



tions. Since x = 1/^—4 Qk (pi[), the physical timescales 
over which the shear evolves are tied to the isotropic cur- 
vature lengthscale, consistent with the previously noted 
interpretation of Type IX anisotropics as maximal wave- 
length gravitational waves. 

Although the temperature anisotropics remain 
quadrupolar, E and B modes mix; in fact PC07 eqs 
(33) or equivalently (50, 52) imply the polarisation 
type oscillates with period Ar] = tt. This effect can be 
traced to the need for the Killing fields to be everywhere 
non- vanishing and yet defined independently of the path 
followed between two points, leading to a spiralling 
of basis vectors perpendicular to a geodesic (see for 
example Fig. 1 in Ref. [38l). 

Overall, type IX models can fit almost any combi- 
nation of quadrupole anomalies in T, E and B modes 
but, since no B modes can be generated while the uni- 
verse remains optically thick, the ratio of the E to B 
anomalous quadrupole amplitudes determines (in the ab- 
sence of reionisation) the conformal time to last scat- 
tering and hence constrains the FRW density parame- 
ters through the Hubble history. Reionisation compli- 
cates this picture, as does the existence of inhomogeneous 
modes superimposed on the FRW background, but in 
principle anomalous quadrupole moments are a powerful 
constraint for observers in closed anisotropic universes. 



V. CONCLUSIONS 

I have discussed anomalous CMB signals arising from 
broken FRW isotropy in open, fiat and closed universes 
as described by mildly anisotropic Bianchi models. These 
are of interest both pedagogically and in the context of 
known CMB anomalies. It is possible to build anisotropic 
models which evade traditional constraints from nucle- 
osynthesis by hiding the anisotropy in super-horizon 
gravitational waves which are pure gauge at early times. 
Whether such modes can be quantitatively linked to 
the observed CMB anomalies remains an open question 
which will be addressed in a future work. 

It is important to stress that all effects discussed in 
this work, including for instance E-B polarisation mix- 
ing, are purely gravitational and require no magnetic 
fields, extensions to general relativity or exotic matter 
source terms. Furthermore the models are predictive - 
polarisation-to-temperature ratios and (in the case of fiat 
and open models) the detailed structure of the possible 
anomalous patterns are dictated by FRW density param- 
eters which are constrained by a range of other observa- 
tions. Thus a future Bianchi signal detection yielding 
consistent cosmological parameters, while highly surpris- 
ing from the standpoint of present infiationary frame- 
works, would have to be taken seriously. 
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